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Dimensionless ratios of physical properties can characterize low-temperature phases in a wide va¬ 
riety of materials. As such, the Wilson ratio (WR), the Kadowaki-Woods ratio and the Wiedemann- 
Franz law capture essential features of Fermi liquids in metals, heavy fermions, etc. Here we prove 
that the phases of many-body interacting multi-component quantum liquids in one dimension (ID) 
can be described by WRs based on the compressibility, susceptibility and specific heat associated 
with each component. These WRs arise due to additivity rules within subsystems reminiscent of 
the rules for multi-resistor networks in series and parallel — a novel and useful characteristic of 
multi-component Tomonaga-Luttinger liquids (TLL) independent of microscopic details of the sys¬ 
tems. Using experimentally realized multi-species cold atomic gases as examples, we prove that the 
Wilson ratios uniquely identify phases of TLL, while providing universal scaling relations at the 
boundaries between phases. Their values within a phase are solely determined by the stiffnesses 
and sound velocities of subsystems and identify the internal degrees of freedom of said phase such 
as its spin-degeneracy. This finding can be directly applied to a wide range of ID many-body 
systems and reveals deep physical insights into recent experimental measurements of the universal 
thermodynamics in ultracold atoms and spins. 


I. INTRODUCTION 

One of the central challenges in condensed matter 
physics is to understand how different phases of mat¬ 
ter can arise and how these phases can be characterized. 
Many phenomena such as, e.g., superconductivity, mag¬ 
netism and quantum phase transitions in strongly cor¬ 
related systems [T], Bose-Einstein condensation of dilute 
gases [5] and of excitons in semiconductors, electronic 
transport in low-dimensional systems and heavy-fermion 
physics [3], are known to exist due to the collective nature 
of the underlying many-body processes. Collective phe¬ 
nomena are particularly strong in low-dimensional sys¬ 
tems where the reduced dimensionality enhances the in¬ 
teraction of elementary constituents gH7j. In order to 
characterize the various phases, dimensionless ratios such 
as the celebrated Wiedemann-Franz (WF) law [5] or the 
Kadowaki-Wood ratio mm are very useful. They usu¬ 
ally involve ratios of measurable quantities which stem 
from similar underlying processes. 

The WF law is universal across a wide range of ma¬ 
terials and temperature regimes because non-universal 
contributions due to, e.g., density of states and effec¬ 
tive mass often cancel out. Conversely, deviations from 
the WF law can be used to characterize the emergence 
of new physical processes nn. Remarkably, recent new 
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experiments )T2j show that the WF law holds even at 
quantum phase transitions. Similarly, the Wilson ratio 
(WR) 0i, 

K = i(^-)‘x (i) 

^ \ M B ^Lande / 


between the susceptibility \ and the specific heat cy di¬ 
vided by the temperature T, measures the strength in 
magnetic fluctuations versus thermal fluctuations m eh 
nn. Here kg is Boltzmann’s constant, fig is the Bohr 
magneton and ffLande is the Lande factor. This dimen¬ 
sionless ratio has been observed in a wide variety of 
Kondo systems 0 HU- Recent studies of the WR for 
magnetic states of a ID spin ladder compound m and 
the two-component attractive Fermi gas [18] show that 
Ryj allows a convenient identification of magnetic phases. 
Dimensionless ratios therefore provide an elegant exper¬ 
imental and theoretical approach to a qualitative and 
quantitative characterization of the nature of complex 
multi-component quantum liquids. 

In Fig. [T] we show that the Wilson-like ratio 
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relating particle fluctuations to energy fluctuations, is 
even more successful in determining phases in quantum 
liquids. Recent experiment [T9] show that the compress¬ 
ibility WR ([2]) determines the Luttinger parameter for 
the phase of TLL in ID Bose gases and characterizes the 
quantum fluctuations at quantum criticality. The key ob- 
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FIG. 1: Phase diagram for the attractive SU(2) Yang-Gaudin 
model at T = 0 . 00 l£ 2 /fcs given as 3D plot of the WR Rfy 
in the p — H plane. The pair binding energy is denoted by 
£2 = 0.5 and the interaction strength <710 = —2. The expected 
values of Rfy = 1,4 in the excess fermion phase (F) and the 
pure pair phase (P) in strong coupling limit, respectively, are 
indicated on the legend on the right. The mixed phase of 
fermions and pairs is denoted as F+P. 


servation underlying the predictive strength of both WRs 
is that the phases are related to simple additivity rules of 
the underlying elementary excitations. In this paper, we 
shall prove that such additivity rules are in fact general 
for multi-component quantum liquids in ID — and we 
expect that they also hold for multicomponent Fermi liq¬ 
uids in higher dimensions m■ We further show that both 
WRs stem from these exact additivity rules, can quanti¬ 
tatively identify quantum phases of Tomonaga-Luttinger 
liquids (TLL) independent of microscopic details and ex¬ 
hibit universal scaling behavior at the quantum transi¬ 
tions between phases. Their values within a given phase 
give information about the internal degrees of freedom 
of a phase, e.g., its spin-degeneracy. Such information is 
important, e.g., for the collective nature of multicompo¬ 
nent interacting quantum liquids in cold atoms l20r l23| . 
large-symmetry fermionic systems [241 W>\ . spin chains 
[T71121)] and ID Hubbard model [5] . 


II. THEORY 


A. General considerations 

In order to show the versatility of the WRs in iden¬ 
tifying phases, we are interested in systems that have a 
rich phase diagram. We therefore start by studying first 
an attractive quantum liquid that supports a hierarchy 
of bound states. The system shall consist of unbound 
particles, pairs of particles, triples of particles and so 
on until we have at most w-tuples. For convenience, we 
shall denote a bound state formed from r particles as 


a r-complex. Let N r denote the number of r-complexes 
which have formed. Then the total number of particles 
is given as N = Ni + ... + wN w and the density of a 
r-complex is n r = N r / L such that the total particle den¬ 
sity is n = rn r- The relative propensity of a phase 

is governed by a set of external fields H r , r = 1,..., w 
with H w = 0. Here the H r are coupling to the (spin) 
moment of each r-complex. Then we can write for the 
Hamiltonian of the system 


W— 1 

H = T + V - J2 H r N r - l lN : (3) 

r= 1 


where T and V are as of yet unspecified kinetic and many- 
body interaction energies. 

In order to compute the WRs ([!]) and ([2|, we need 
to compute the susceptibility y, the compressibility n 
and specific heat cy for the urcomponent system in the 
Gibbs ensemble G(p, Hi,..., H w ). In a single compnent 
system, y and n can be straightforwardly computed as 
dM/dH (with magnetization M ) for constant particle 
number (canonical) or dn/dp for constant external fields 
(grand canonical). For a multi-component system, we 
therefore define in complete analogy a chemical potential 
p r for each of the r-complexes via 

1 £f 

Hr = fX - 1 — H r -\ -, (4) 

r r 

where e r denotes the binding energy of an r-complex. 
Remarkably, the quantity P 2 has already been measured 
in a recent experimental study of the equation of state for 
2D ultracold fermions [23]. It gives a deep physical in¬ 
sight into the crossover from Bose-Einstein condensate to 
Bardeen-Cooper-Schrieffer superconductor. Indeed, the 
choice Q allows us to define a Fermi energy at T = 0 
for each fluid of r-complexes in the same way that is 
defined in the Landau’s Fermi liquid picture. We now 
introduce a stiffness in grand canonical and canonical 
ensembles, respectively, as 


d: 
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for a r-complex fluid subject to /1 and the field H r t. Also, 
the sound velocity for r-complexes is defined as usual via 
v r = | k F , where e r (fc) is the dispersion relation for 

r-complexes and kF the Fermi momentum. With these 
definitions, we can then prove in ID that the individual 
K r , xr and cy, r in terms of the Df, D x , and v r are given 
as 


r 2 1 r 2 1 7 ik 2 B T 1 

7r hDf’ ^ r,r 7 T%D X y' ° V,r 3 H v r 


as shown in the Appendix. We note that similar relations 
hold in Fermi liquids (FLs) [T6|. Our strategy therefore 
is to derive the densities n r , the chemical potentials p r 
and the dispersions e r (k) for the r-complexes as functions 
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of Hi ,..., H w and g. We expect that the definitions ^ 
and ([h]) are useful in general for interacting quantum liq¬ 
uids. However, analytical or numerical access to these 
quantities is not necessarily straightforward. Following a 
series of recent papers m , we prove here that integrable 
ID multi-component systems allow the explicit construc¬ 
tion of n r , g r and e r using the thermodynamic Bethe 
Ansatz (TBA) Following earlier theoretical [Hj 

and recent experimental mm 2 findings, we furthermore 
speculate that similar relations also hold for FLs in higher 
dimensions and expect that such effective chemical po¬ 
tentials for multi-component systems can therefore serve 
as convenient handles to describe multi-component quan¬ 
tum liquids. 


B. The w-component quantum liquid 


Let us now consider a specific w-component Hamilto¬ 
nian ([3]) that supports multi-component quantum liquids. 
A convenient and experimentally relevant example is the 
ID SU(w) Fermi gas with (5-function interaction confined 
to length L Ep S3. The system is described by the 
Hamiltonian ((3J) with 


T+V 


2+fflD £ Sfr-xj), (7) 

i =1 1 1 <i<j<N 


and with the chemical potential /i and the effective Zee- 
man energy E z = ]C”=i \r{w—r)n r H r . There are w pos¬ 
sible hyperfme states 11 ), | 2 ),..., |w) that the fermions 
can occupy. Experimentally, gm (= — 2fi 2 /maiD, with 
aiD the effective scattering length in ID [33]; from now 
on, we choose our units such that tr = 2m = 1 un¬ 
less we particularly use the units.) can be tuned from a 
weak interaction to a strong coupling regime via Feshbach 
resonances. This model provides an ideal experimen¬ 
tal testing ground to probe few- and many-body physics 

E3EZHM1- 

For the system described by ([7| , the relation Q follows 
naturally from the structure of the TBA equations [27] 
(see Appendix[A|. We will use in addition general energy- 
transfer relations for breaking a w-complex into smaller 
r-complexes, i.e. 

— H r = fl r g w T £ w ~£r (8) 

r w r 

with r = 1, 2,..., w — 1. Note that e r = ^r (r 2 — l) g 2 D 
is the explicit binding energy for each r-complex in the 
system 0 -, Using the energy and particle conservation 
conditions, 0 and n = i rn r, respectively, we find 
that k and \ r obey the additivity rules, 


K 
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Xr 
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in the TLL phases (see Appendix H- Here the suscep¬ 
tibility Xr represents the responses of different bound 


states to the change of the field H r . Such additivity ap¬ 
pears naturally for a non-interacting fluid. Determining 
v r from the TBA, we similarly find that 


c v = c v,i + c v,2 + • ■ ■ + Cy )V 


( 11 ) 


Based on the rules given in 0 . ( fl 0 |) and flTTj ), we can 
construct the WRs of the w-component SU (wjsystem to 
be 
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( 12 ) 

(13) 


These ratios are dimensionless and uniquely determined 
by the sound velocities and stiffnesses. We note that the 
form of the WRs in (12) and (131 is similar for a w- 
component FL; all interaction effects have been included 
into ([ 8 ]) via the choice of g r . 

In the strong coupling regime, for a pure r-complex 
phase (i?^r = 0 in pure phases fTS]) can be given in the 
form 


R^ r = vK r = r 2 (l - 2 \B r ±- + B 2 ^j + 0( 7 " 3 ),(14) 

where 7 = giy,/2n is the dimensionless interaction 

r— 1 

strength, B\ = 0 and B r = 1/fc with r = 2 ,..., w. 

k =1 

In the above equation, I\ r is the phenomenological Lut- 
tinger parameter which can be directly measured through 
the Wilson ratio R This provides an important way 
to test the low energy Luttinger theory. For 7 — 00 , Rfy 
then displays plateaus at the integers l 2 ,..., r 2 ,..., w 2 . 
Deviations from integer r 2 = 1,4,9, occur not because of 
a weakness of the ratios introduced here, but because of 
the physics underlying and determining the stability and 
mixture of bound states in a given phase. Thus Rfy pro¬ 
vides a convenient quantitative phase characteristic for 
quantum systems at finite T as well as at T = 0, see 
Fig □ and below. 


C. Scaling of the WRs 


It is particularly interesting that the WRs identify non- 
TLL behaviour and quantum criticality in the quantum 
critical regime. In fact, the WRs show sudden enhance¬ 
ment near a quantum phase transition due to a break¬ 
down of the quantum liquid nature, i.e. the vanishing lin¬ 
ear dispersion for ID systems. For both WRs, we have 
the scaling law 


p«,x _ 7r«,x 


(V-Vc) 
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iv - Vc) 
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,( 15 ) 


where j3 = ( d/z ) + 1 — 2,/{vz) with z = 2, v = 1/2 
and d = 1 for ID systems is universal; F K ' X , Q K,X are 
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III. APPLICATIONS 


A. The Yang-Gaudin model. 


In order to show the usefulness of the new ratios, we 
start by considering the w = 2 ID spin-1/2 Fermi gas 
with a ^-function interaction [.'151I3S] . This Yang-Gaudin 
model has already provided an ideal experimental test¬ 
ing ground [33133133 for few- and many-body physics 
[23 [391140]. For Ryy, the additivity rule of susceptibility 
and its connection to the 74/’s following ( 12 ) have been 
reported in [IS]. For 77/y, we now find from (13) that 


FIG. 2: Scaling of Ryy at phase boundaries (a) V-F and 
(b) V-P with H/t 2 = 1.16 and 0.9, respectively (cp. Fig. [1]. 
The left panels in (a) and (b) show Rw vs. g/t 2 for reduced 
temperatures t = 0.003 (A), 0.002 (□) and 0.001 (o). The 
vertical dotted lines indicate /u c /e 2 . The dashed lines show 
the calculated slopes C r /t. The right panels show the scaled 
Rw as in ( |15| ) with /? = —1/2 vs. (/.i — g c )/t 2 1. The vertical 
dotted lines are as in the left panel while the horizontal lines 
indicate the strong interaction limit of Ryr in the (a) F and 
(b) P phases. 


the scaling functions (see Appendix [C]). The second term 
in © with /3 = —1/2 reflects a contribution from the 
background at the temperatures above the energy gap 
A ~ \r] — r] c \ zv , with driving parameter 77 (e.g. g, or H r ). 
These critical exponents agree with the experimental re¬ 
sult for the thermal and magnetic properties of the ID 
spin chain at criticality (19; 20j ■ For a phase transition 
from the vacuum to a TLL phase, Ao = 0. Based on 
(15), the slope of the temperature-rescaled Ryy curves at 

a critical point is given by ^g^ ^J = C r /T , see Ap¬ 

pendix C. 



(16) 


Writing the explicit dependence on fi, H and g id is pos¬ 
sible, but tedious (we give explicit forms for 74/ 2 and 
Vi t 2 in Appendix E- In Fig. |T] we show R^ in the g,- 
H plane. Ry, elegantly maps out the three finite tem¬ 
perature TLL phases of 1-conrplexes, i.e. fully-polarized 
fermions (F), 2-complexes, i.e. pairs (P), and the mixed 
TLL of excess fermions and pairs (F+P). The empty vac¬ 
uum phase (V) is also present in Fig. |TJ In particular, 
we note that for the pure phases F and P, R^y ~ r 2 , 
i.e. close to the strong coupling value rK r m■ Hence 
for the Yang-Gaudin model in Fig. |T]we see R^ = 1 in 
phase F and 4 in phase P. In mixed phases, e.g. F+P, 
such a constant plateau no longer exists. For R^y in the 
critical regime, we find a rapid increase due to a strong 
increase in thermal fluctuations. In Fig. [2] we show the 
scaling behavior (15) for 7?/y close to the transition from 
the vacuum phase into the F and the P phases. We find 
Ci « 1.2546 and C2 ~ 2 x 5.0185, respectively. Eq. 
(15) does not contain any free fitting parameter, so one 
can use this scaling to determine the temperature of the 
quantum liquid. The scaling law is also similar to the 
scaling in the non-FL regime of heavy fermions 131 132 ] • 


The significance of the WRs can also be under¬ 
stood from the quantum fluctuations of magnetization, 
( SM 2 ) = £ d kBTx, and particle number, (5N 2 ) = 

^ksTn ; l d denotes the observation volume in d- 
dimensions. Therefore, microscopically, y and k measure 
the strength of these fluctuations just as cy quantifies the 
energy fluctuations; macroscopically, the temperature- 
independent compressibility and the linear-temperature- 
dependent specific heat remarkably preserve the nature 
of quantum liquid at the renormalization fixed point, see 
Eq. (16]) like the Fermi liquid m- Consequently, 7?^ and 
7?// characterise the competition between fluctuations of 
different origin. A constant WR implies that the two 
types of fluctuations are on an equal footing, regardless 
of the microscopic details of the underlying many-body 
system. On the other hand, the growth of the WRs in the 
critical regime indicates the long-range character of the 
quantum fluctuations at the quantum phase transitions. 


B. A w >3 quantum liquid. 

An even richer phase diagram exists for the three- 
component 5-function interacting Fermi gas with an at¬ 
tractive interaction m- Its T = 0 phases are known 
to consist of excess fermions (F), pairs (P), trions (T) 
and mixtures thereof H? ]T. Using the TBA equations, 
we again numerically and analytically calculate the WRs 
and determine their Hi, H 2 and T dependencies. In 
Fig-13 we show that in the pure phases F, P and T, we 
have R 1 ^, = 1,4,9, respectively, and the phase bound¬ 
aries are clearly marked by large increases in Ry? near 
the critical points. At T > 0, we see that the n and 
cy curves become progressively more rounded across the 
phase transitions for T increases. The magnetic field as¬ 
sociated with the position of the finite-height peaks varies 
as a function of temperature. This can be used to define 
a crossover temperature T* for each such magnetic field 
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FIG. 3: Upper panel: plot of Rw for the ID three-component 
attractive Fermi gas in the Hi—Hi plane at temperature t = 
0.001E3 /ks- The values of Ryy = 1,4,9 in the pure phases 
F, P and T, respectively, are indicated by the colors marked 
on the color scale shown. The dashed lines follow the phase 
boundaries as indicated. Lower panel: compressibility k and 
specific heat cv IT vs. external field Hi I £ 2 , for a fixed choice 
of polarization m = ni indicated by the dash-dotted line in 
the phase diagram (upper panel). At zero temperature, n and 
cv/T satisfy the additivity rules (j9j and (111 as shown by the 
solid lines while at finite T, they exhibit peaks as indicated 
by the arrows for t = 0.003e3/fes. Near the phase boundaries 
both k and cy /T diverge quickly as (Hi — IT^) -1 ' 72 . 


corresponding to a peak. For T < T* , our analytical and 
numerical results reconfirm the validity of the additiv¬ 
ity rules © © i-e- k and x remain independent of the 
temperature and cy depends linearly on the temperature 
in the TLL phase (cp. lower panel in Fig. [3]). This nature 
is also seen in Fermi liquids mm]. The additivity na¬ 
ture of the susceptibility is presented in Appendix [B] The 
significance of such additivity rules is the characteristic 
of quantum liquids at the renormalization fixed point, 
such as TLLs in integrable models and Fermi liquids in 
metals, etc. For T > T* , Rtfo exhibits the critical scaling 
behaviour (15) (see also Appendix [C|). Fig. [4] shows the 
contour plot of Ryy at low temperatures. The univer¬ 
sal quantum critical behaviour © is characterized by 
the exponents z = 2, v = 1/2. For T < T* , the critical 
exponents are z = v = 1 . 


Moreover, using the TBA equations for the SU(w) 
Fermi gases with an attractive interaction [45] , we calcu¬ 
late the susceptibility \r — dM/dH r in response to the 
change of magnetic field H r for the spin-gapped, phases. 
The other magnetic fields are fixed and the magnetization 
is given by M = X^li 1 n r r(w — r)/2. In dimensionless 
units, the energies are rescaled by the interaction energy 
eb = \h 2 c 2 /(2m), i.e. the temperature t = T/e&, the 
pressures p r = p r /(\c\eb), the susceptibility x = |c|£/ 2 , 
see the Appendix B. We obtain explicitly a general ex¬ 
pression of the susceptibility Xr for the gapped phase in 
which a small number of the r—complexes are created 



FIG. 4: Mapping out quantum criticality with R^ in the T— 
Hi plane. The Hi values have been chosen to follow the fixed 
polarization m = ni as indicated by the dash-dotted white 
line in Fig. [3] The white dashed lines denote the crossover 
temperatures T* beyond which the TLL phases break down. 
The inset shows individual temperature curves of Rw vs. Hi 
for representative temperatures. The bold solid line shows 
the result corresponding to (131 while the thin dashed lines 
show the numerical results obtained from the TBA equations 
(Appendix [A]). 


due to the change of H r for t -C 1, 

1 1 


Xr 


4\[2/k \Jt 

where the gap is given by 


yjr(w — r)re * 


w min(r,ra) 


A. r « —rp r + ^2 ^2 


4p n 


m—1 q 

2 q ^ r + 


' m(r + m — 2 q)' 


(17) 


(18) 


see the Appendix C. For example, from the upper panel 
phase diagram in Fig. [3j we observe that in the gapped 
phase of T, xi with Ai = — p,\ +2p 3 /3 and X 2 with A 2 = 
2 / 7,2 + 16p 3 /9 show a dilute magnon behaviour related 
to the phase transitions from the phase T into the mixed 
phase F+T and into the the mixed P+T, respectively. 
They decay exponentially as the temperature approaches 
zero. By properly choosing H\ and # 2 , we also can have a 
phase transition from the gapped phase T into the mixed 
phase F+P+T. 


C. Critical theory for the SU(w) repulsive Fermi 
gases 

The existence of internal degrees of freedom signifi¬ 
cantly changes the quantum magnetism and the dynam¬ 
ics of the system compared to the spinless Bose gas. 
It is well established that the critical behaviour of the 
spin SU ( w ) chain can be described by the Wess-Zumino- 
Witten a model of level £ = 1 with the Kac-Moody cen¬ 
tral charge C s = £(w 2 — 1)/ (£+w) [49]. The w-component 
repulsive Fermi gases display a U(l) <S> SU(w) symmetry 
characterized by one charge degree of freedom and w — 1 
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spin rapidities. The low-energy physics of the system 
is described by the spin-charge separated conformal field 
theories of an effective Tomonaga-Luttinger liquid and an 
antiferromagnetic SU(w) Heisenberg spin chain [51I5U]. 
By using the TBA equations given in [35; with H —> 0 
we find that the pressure of the ID repulsive Fermi gases 
is given by 


P = Pt =o 


w(w 2 — 1)H 2 
247 TVs 


nT 2 r 1 C s - 
6 Vv c 


(19) 


where pr=o is the pressure at T = 0 and v a , c are the 
pseudo Fermi velocities in the spin and charge sectors, 
respectively. This result is consistent with the critical 
field theory for the SU(w) spin chains QHj. In the zero 
magnetic field limit, the spin and charge velocities for the 
Fermi gas with strong repulsion are, respectively, 


2n 2 7r 3 aiD n , Q v ^ 

v s = --- (1 + 3ai D l? n). 

ow 

v c = 2 n 7 r (1 + 2amZn) ? 


( 20 ) 


where Z = [i’iw) + C] and C the Euler q-constant, 
ip(x) the Euler ^-function. The susceptibility is then 
given by Luttinger-liquid relation 


XVs = w(u> 2 - 1 ) 


( 21 ) 


in the limit H —> 0. Indeed the sound velocity v c (20) 
of the Fermi gas in the large-w limit coincides with that 
for the spinless Bose gas, whereas v s vanishes quickly as 
w grows. This gives a reason why the quantum liquid 
of multi-component fermions reduces to the liquid of a 
spinless Bose gas in this limit [3J. In this sense, Ryy 
captures this unique large-spin charge separation mecha¬ 
nism in the w-component repulsive Fermi gas. Its explicit 
expressions is 


w(w 2 — l)v c 
3 [(w - l)v c + v s ] ’ 


( 22 ) 


displaying plateaus of height w(u> + l)/3 for either strong 
repulsion or in the large-w limit, hence capturing the spin 
degeneracy. For example, Ryy = 2 and 4 for the two- 
and three-component Fermi gases with strong repulsion, 
respectively. 


IV. EXPERIMENTAL REALIZATIONS 

Both WRs Ryy and Ryy are readily accessible by ex¬ 
periments. For example, the ID SU(u;) (5-function in¬ 
teracting Fermi gas of 173 Yb has been realized experi¬ 
mentally [M] . It was shown that in the large-w limit 
the ground state of the gas with a repulsive interac¬ 
tion exhibits properties of a bosonic spinless liquid. In 
the context of ultracold atoms, it is highly desirable to 
measure the quantum criticality and the TLL in such 
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FIG. 5: Upper panel: Contour plot of the WR R%y Q for 
the Yang-Gaudin model in a harmonic trap with fixed par¬ 
ticle number N and three different polarizations P = n\/n 
at the temperature T = 0.00l£2 /ks- The WR evolves thin 
round pe aks near the phase boundaries, whereas the values of 
the WR ( |16| ) quantify different quantum phases of the TLLs. 
The WR peaks are in good agreement with the experimen¬ 
tal phase boundaries (red and yellow dots) observed in [37| . 
Lower panel: the red-dashed and yellow-solid lines show the 
WR Ryy in the trapped gas with the polarizations P = 0.027 
and P = 0.341, respectively. They indicate segments of the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in the centre 
accompanied by the wings of the P and the F states, respec¬ 
tively. The highest Ryy plateau approaches 4 suggesting a 
quasi-ID superfluid nature, whereas the lowest plateau shows 
the free fermion nature with R^y = 1. For a critical polariza¬ 
tion P c = 0.176, the green-dashed line shows that the trapped 
gas consists solely of the FFLO-like state. 


Fermi gases with rich spin and charge degrees of free¬ 
dom. In this scenario, the Yang-Gaudin model is an 
ideal model to conceive critical phenomena induced from 
spin and charge interaction effects. This model was 
recently studied via an ultracold atomic gas in a har¬ 
monic trap, such as the two-component ultracold 6 Li 
atoms of Ref. Due to the harmonic confinement, 

the chemical potential in the equation of state should 
be replaced by p(x) = p ( 0 ) — muj 2 x 2 /2 (within the lo¬ 
cal density approximation). Here x denotes the posi¬ 
tion along the ID trap. Changing x is then equivalent 
to changing p and different phases are located at differ¬ 
ent spatial positions along the trap. Using a rescaled 
coordinate y = xa\Y>/(2a 2 ) with axial characteristic os¬ 
cillation length a = \Jh/mu) , the density profile of the 
trapped gas can be determined from the dimensionless 
quantities Naf D /a 2 = 2oid /_ dyn(y) and polarization 
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Maf B /a 2 = 2aiD J^° 00 dym z (y). The compressibility can 
be extracted via re = dn/dfj, = —dn/dx/(xmtu 2 ) [ 151117 ] , 
while the density n(x) is read off from experimental pro¬ 
files along x. The specific heat can be measured through 
the sound velocities, while the density waves of pairs and 
unpaired fermions can be experimentally created by a 
density depletion or pulse with a far-detuned laser beam 
[48] . In the upper panel of the Fig. [5j we observe that 
R^r at a very low temperature naturally maps out the 
phase diagram of the system, showing a good agreement 
with the experimental phase boundaries [37] . In the lower 
panel of Fig. [5j we show the behaviour for 1?^, exhibit¬ 
ing the plateau structure as the system passes through 
different phases upon changing x. 


V. CONCLUSIONS 


In summary, we have shown that the additivity rules of 
physical properties reveal an important characteristic of 
the TLL and provide the physical origin of the dimension¬ 
less WRs. Such dimensionless ratios can be used to iden¬ 
tify full TLL phases and capture the essence of phenom¬ 
ena ranging from quantum criticality to spin and charge 
separation in a wide variety of ID many-body systems. 
We have presented some universal relations, such as the 


WRs (12 1 and (13), the Luttinger parameter (14), the di¬ 


mensionless scaling function of the WR ([15]) , the suscep¬ 
tibility for the gapped systems ( fl7| ) and the WR of SU(w) 
repulsive Fermi gas related to the level-1 Wess-Zumino- 
Novikov-Witten conformal theory ( [22] ). We also show ex¬ 
cellent agreement between the phase diagram predicted 
by the WRs and the experimentally determined one for 


the Yang-Gaudin model [37]. Our results therefore suc¬ 
cessfully demonstrate how to predict universal laws for 
experimentally realizable quantum liquids in Bose and 
Fermi degenerate gases, Bose-Fermi mixtures, ID Hub¬ 
bard model, strongly correlated electronic systems, spin 
compounds near to and far from quantum critical points 
on an equal footing. 
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Appendix A: Thermodynamic Bethe Ansatz for the 
SU(ui) Fermi gas 


The Gibbs free energy of the SU(w) attractive gas is 


rT 


G(/b Hi,..., H w ) = ^ — / dk In 1 + e ~ e ^ T 


r—1 


with dispersions e r (k) defined in the TBA via 


(Al) 


r—1 ( w w 

e r (k) = rk 2 — r/i — H r — e r + ^ ^ o g+r _ 2p * Tin 1 + e~ Cq< - k ^ T + ^ a q _ r *T In 1 + e~ eq< - k ^ T 

p= 1 Kq=P q=r+ 1 


— '^2 a q * T In 


9=1 


1 + e 


-Vr, q /T 


ijr :l (k ) = l ■ (2 H r - H r _ 4 - H r+1 ) + ai*T In [l + e~ eAk)/T ] + U iq * Tin |~1 + e -^,«(fc)/r 


9=1 


- s k * nn [l + e ~ riT - 1 - q{k)/T ] - Siq * Tin [l + e ~ rir + 1 ’ q{k)/T 
9=1 9=1 


with 


Uij (x) = 
Sij(x) = 


a\i~j\(x) + 2a\i_j\ +2 (x) + ... + 2a t+j - 2 (x) + a i+j {x), l ^ j 
2a 2 (x) + 2a 4 (x) + ... + 2a 2[ - 2 (x) + a 2i (x), l = j, 

a \i-j\+i(x) + 2a\ l _ j \ +s (x) + ... + 2a i+j _ 3 {x) + a /+J -_ 1 ( 1 ), l ^ j 
ai(x) + a 3 (x) + ... + a 2 i- 3 (x) + a 2 i-i{x), l = j, 


(A2) 


(A3) 


r 


while * denotes the convolution (a * b)(x) = J a(x — eters; furthermore a n {x) = n|c|/27r[(nc/2) 2 + x 2 ]. Here 
y)b{y)dy and the r) r j{k) represent the spin string param- 























c = mgm/h 2 = — 2 / oid - The numerical results used 
in the figures were obtained by solving the above TBA 
equations. 


Appendix B: Additivity rules for SU(w) Fermi gases 


From Q, we find that constant fields in the grand 
canonical ensemble, i.e. dH i = dH 2 = ... = dH w -\ = 0, 
imply dg r = dg for r = 1,..., w. Consequently, we have 


r\ \—\W W r\ W 

an ron r ^ 

= dii = dg = ^ ^7 = 2^ K ' 

r —1 7"—l 


(Bl) 


with n r as given by ([5| and ([6]). Note that this result is 
general and does not use special properties of the SU (w) 
gases. 

The additivity rules for the spin susceptibility are very 
intriguing. For our convenience in analysis of the SU(w) 
case, we prefer to keep the ratios among the densities of 
the charged bound states as ni : n 2 : ■ ■ • : n w - 1 = 
Ai : A 2 : ... : A^-i, then we can parameterize the 

^1 7 ^2 ? • • • J 

n r = -Tfi, (r = 1, 2,..., w — 1), n w = — (n — n) (B2) 
A w 


explicitly: H r = D r E and the susceptibility in response 
to the H r can be expressed as 


1 

A> 


d H r 
d M 


YXl K-Dr 




• T> r E. 


We note that the last term in the Eq. (B3) linearly de¬ 
pends on the density n r and n w . Therefore it can be 
safely dropped off in the calculation of the susceptibili¬ 
ties because of the second order of derivatives. In terms 
of the ground state energies E = EE l f° r indi¬ 
vidual charge bound states We define the stiffness as 


1 

Xrt 


= —s-D™, = 


Er=l *rV r 


EEi \r{w-T)X r 


■ D r E t , (B 6 ) 


where r = 1,2,..., w — 1 and i = 1, 2,..., w. Conse¬ 


quently, we have the susceptibility ( 10 ) in responses to 
the external field H r . 

In fact, Zeeman spliting can be characterized by the 
Zeeman energy levels e r z or by the effective magnetic 
fields H r with r = 1,2, Here H w = 0. Both 

sets of parameters are related via the relation 


= — H r (n r — n r+1 ). 


(B7) 


where A = Ai + 2A 2 + ... + (w — l)A u ,_i. In order to 
compute the additivity rules for x ” 1 = dH r /dM , with 
magnetization M = Y^r=i r ( w ~ r ) n r/ 2, we start with 
the Legendre transformation from Gibbs free energy to 
ground state energy: E = G + fm + EEi* n rH r , and 
in the ensemble {ni,--- ,n w -i,n} the field H r can be 
obtained: 


Hr = 


dE 


dn, 


dE 

dn r 


n i,--- ,n r _i,n r _|_i,••• ,n w — i,n 

r dE 

ni ,••• ,n r — i ,n r +i ,••• ,n w W dn r 


ni,-" ,nw — i 


(B3) 


where E = EEi ^ r denotes the energy of the multi- 
component ground state m and we have dE/dn r = rg r 
for consistency with ([4]). We emphasize that the additiv¬ 
ity of E is a fundamental property of a TLL, implied by 
the linearity of the dispersions [ 171 - 
Following the relation (B2 1 , we can define differential 
forms of any thermodynamic function / = /(m, • • ■ , n w ) 


, f V df Xr A ~ ° f 1 

d/ = 2 ^ «-v-An “ n- 

f on r A on w w 

r—l 

u, - 1 i \ 

AM = ^-r(w — r)-T • An. 

r—l 


An, 


(B4) 

(B5) 


We denote the operator V r = with 

L ' OTlr W OTl w 

r = 1,2,..., w — 1. Then the field H r can be expressed 


A consistent solution of this equation gives the relations 
between H r and Zeeman energy levels [45j . If we denote 
the difference between the energy levels of fermions in 
the states |m + 1 ) and |m) as A m+ i im = e ^ +1 — e™, 
then the total susceptibility x = dAt\ i wdd Atotal = 
EEa A r+ i ir is given by 


X 


-l 


dE 


w —1 
r—l 


A 


r+l,r 


AM 

1 1 

Xi Xw-i' 


d {H i + H w _ i) 

AM 

(B 8 ) 


Using the numerical calculation for the Gibbs free energy 
G(n, Hi,H 2 ) from Eq. ( |Al[ ), we observe that the suscep¬ 
tibilities Xi an d X 2 satisfy the relation (10 1 in the TLL 
phase of trions, pairs and single atoms, see Fig. 

Using the TBA equations [251 US] , one can easily prove 
that the specific heat can be written in terms of sound 
velocities, i.e. 


Cy 


7 T 2 fc|T y. 1 

3 fvKVr 

r 


(B9) 


Appendix C: Scaling forms on boundaries 

We focues on the low temperature behaviour for w = 3 
to derive the analytical results for the thermodynamics, 
and finally we derive formula of general SU(w) gases of 
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X\ X 2 



FIG. Al: Susceptibilities in response to H± and H2 for the three- 
component Fermi gas with a fixed ratio of ni/712 = 1 at finite 
temperatures. Under this setting the H\ and if2 have a one-to- 
one mapping along the line of ni/712 = 1 with a fixed total den¬ 
sity. These two figures show that the susceptibility is temperature- 
independent for the TLL phase of F+P+T indicated by the long 
dashed line in Fig. [3] in the main text. The deviation from the solid 
lines in the two figures show a breakdown of the TLL. Here T, B 
and F stand for trions, pairs and excess fer mion s, respectively. The 
solid lines confirm the additivity rule Eq. {To} in the main text. 


the cricarity about susceptibilities.In the strong interac¬ 
tion limit, we can solve these TBA equations analytically 
at low temperatures T £ 3/^3 by simplifying 


e r (k) 


rk 2 — A r . 


r = 1,2, 3, 


(Cl) 


with 


Al = " + Hl + 

+ ^ji y 3,| + Te-W'-^VTe-^IM/T) 


1 9 4 

A 2 — 2/1 + ~c H- H 2 — 


1 16 

c\ Pl ~ \c\ P2 ~ ^c\ P3 


+ |C|3^ 


1 


7 Y 1 


224 


r Vo 


4|c | 3 2 ’5 243|c| 3 3 d 

+ Te- {2H2 ~ Hl)/T e~ j2/T I 0 {J 2 /T), 

, n „ 2 2 8 1 

A 3 = 3/1 +2c - -~pi - —-p 2 - -r-rP 3 
|c| 3|c| |c] 

1 28 1 

^a -)- — —Y 2 a -f- -777—77V 


2 c 


27|i 


•I3 J 2,7 


16|c 


3^3,1 


(C2) 


with Y r>a = —\fjp;T a Li a (—e Ar/,T ). Here the effec¬ 
tive spin coupling constant is J r = 2p r /r\c\ with c = 
mg 1D /h 2 . The polylogarithm function is defined as 

Lin(z) = Er= 1 W and W = Er=o (5TF ( x /2) 2k ■ The 
effective pressures p ri with r = 1,2, 3, of excess fermions, 
pairs and trions, respectively, can be expressed as 


pi = L.f 

p 2 = Y 2 3 
P3 = E 3 


ttV, s 
3 5 


3|« 


iis y 3. 


T 3 L 3 + 

■3 1 ’ 2 


4|, 


V 9 3 + 


112 


•13 


. Yo 


3| 


1 E.s + + 


81|c| 3 3 ’5 
1 


■13^1,* 


81c 


8|c 


r To 


1. Vacuum — Pair 


The above simplified TBA equations (Cl) can be used 
to derive universal low temper properties of the SU(2) 
and SU(3) Fermi gases. In the following, we derive the 
scaling forms of the Wilson ratios in the critical regions 
for the two-component Fermi gas. Here we will use the 
dimensionless units as being explained in the main text. 
The phase boundary for the phase transition from V to P 
phase are p c = — \ for h = H/e 2 < 1- Near this critical 
point, the scaling forms of specific heat, compressibility 
and susceptibility in dimensionless units are given by 



2 gl 

' 2 (A- 

Pc)' 

\c\t 


. t 


k ~ 

- >1 

(2 (M- 

Ac) 

■\/ 7T t 

{ t 



X ~ 0 . 


(C3) 

(C4) 

(C5) 


where the functions Q{x) = ^Lia (— e x ) — \xLii (— e x ) + 
|x 2 Li_ 1 (— e x ) and 7(x) = Li_i(— e x ). 


2. Vacuum — Fully Polarized 


The phase boundary for the phase transition from V 
to F phase are p, c = — |, for h > 1; the scaling forms in 
dimensionless units are 


cy 

\c\t 

k 


X 


-j=n (^=M, 

2V2 Vt V t J 

1 T ( A - fic \ 

2 V t J ’ 

1 T ( P~Pc \ 

8 ^ 2 Vt \ t ) ’ 


(C 6 ) 

(C7) 

(C 8 ) 


where H(x) = |Li| (— e x ) — xLi 1 (— e x ) + x 2 Li_ 1 (— e x ). 


3. Pair — Partially Polarized 

The critical fields for the phase transition from P to 
F+P phase are fl c = —| + ^(1 — h)?, for h < 1; the 
scaling forms are 


cy_ 

\c\t 

k 


X 



Ko1 uWt 



A - A<A 

t J’ 


7 


f At ~ A 
V t 


A - Pc\ 

t J’ 



(C9) 

(CIO) 

(Cll) 


where b = (1 - h) (l + f JY=h) , A, = 1 + ^ 

A 2 = 1-^ + fl, Koi = ^Ai and K(x) = |Li|(-e x )- 

xLii (— e x ) + x 2 Li_i (— e x ). 

2 2 
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4. Fully — Partially Polarized 


For the phase transition from F to F+P phase, we have 
p c = — \ + ^(h — l) 5 , h > 1; with scaling forms 


cy 

\c\t 


1 

2 v / 2 t rt 
Ko 3 - 

Xo4 - 


S 


^3 


2(A Me) 
t 

(2{p-p c ) 


--T 


2\f2ixi V 

A 4 


t 




( 2 (/x - p c ) 


8-\/27 rt V £ 


(C12) 

(C13) 

(C14) 


where the constants are given by a = (h — 

i)(i + &V£=i), a 3 = 4^(1 + ^^^), a 4 = 

8V2a _ 1 and v - 1 ( 1 _ VIA 

W 2 , k o3 - 2v ^ w%/s , and Xo 4 - 7r )■ 

The dimensionless function is given by S( x) = 

A^Lig (—e x ) - v^a^Li 1 (—e x ) + y/2x 2 Li_ i (-e x ). The 

slopes of the Wilson ratio curves at the critical point 

p = p c reveal a unique temperature-dependent feature, 

( dR K \ 

~d^~ ) 

given as 


= Tp is 


Cr 1 - 


ttt 2 (J-'(O)^(O)-J-(0)^'(0)) 

Q( 0) 2 


(C15) 


i.e. is a constant for the phase transition from vacuum 
into an r-complex TLL phase. 


5. Susceptibilities for the gaped phase in the 
attractive SU(w) gases 

The total polarization of SU(w) gas can be expressed 
as rh = £“ = i \n r r(w — r), then we can obtain the sus¬ 
ceptibility: 


- dm ^ 1 d 2 p^ 

y r = —= > -k(w — k)- 


dh k dh r 


(C16) 


in the limit of t —> 0 and |c| >> 1, the leading behavior 
of the second order derivatives of pressures reads 


(C17) 

dh k dh r 2 2 ' 1,k ’ v ' 


Substituting (C17) into (C16), we arrive at the explicit 
form of the susceptibility Xr corresponding to the field 
H r when the other fields are fixed 


t 2 r , , T . , aW 

Xr ~-^^Vr(^-Hr Li _i(-e ‘ ) 

« —=sjr[w -r)re * , 

4y 2n 


(CIS) 


where the gap A r is related to the effective chemical po¬ 
tential through A,. = — A^ r \ which can be determined 
from the TBA equations 


uj min(r,m) 

i (r) = rp - Y, 

m= 1 q — 1 

2 q ^ r + m 


4 p m 

m(r + m — 2 q) 


(C19) 


Here the second approximately equal in (C18) holds when 
A r > 0 which implies the existence of the gap and the 
susceptibility presents exponential decay, otherwise Xr ~ 
4 CbT/=^ = 1 which is a positive constant when t —> 0. 
Note that the result in (C18) is different by a factor 4 
due to the convention of Hi = H/2 in the SU(2) case. 


Appendix D: Explicit formulas for the SU(2) Fermi 

gas 


In the mixed phase for SU(2) Fermi gas with one 
magnetic field H\, we have n = 2n 2 + «i and xi,i = 
(/4 B fl) 2 {dm/dpi) n , Xi ,2 = 2(^ B fl) 2 {dn 2 /dp 2 ) n , with the 
stiffnesses in canonical ensemble D* = A- {dpi/dni) n , 
D 2 = ( dp 2 /dn 2 ) n ■ Similarly, «i = (dm/dpi) H 

and k 2 = 2 {dn 2 /dp 2 ) H with D% = A (dp 2 /dn 2 ) H , 
D* = Y (dpi/dni) H defined in the grand canonical en¬ 
semble. 

For strong coupling, we find the explicit form of the 
p r in terms of the densities of pairs n 2 and unpaired 
fermions ni (in units of h 2 /(2m)) 05] to be 


M2 


Mi 


2 Tin 


. , 2 . njm 

4 + 3|c| + |c| 


+ 


4n\ii2 16 n\n 2 




2 rd 

3M 


4 nf 2>ri\n\ bn\ 
3|c| c 2 4c 2 

(Dl) 


48 n\n 2 4n\n\ 


(D2) 


Here c = mgm/h 2 = — 2/aip is the interaction strength. 
Therefore, the compressibilities and susceptibilities are 
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given by (in units of h 2 /( 2 m)), e.g. 


,-i 


-i 


-l 


7 T 2 n 2 


1 + 


6ni 4n 2 


24 ?r 2 


2 c Til 


24n 1 n2 17n| 2n.2 


27r 2 ni 1 + 


12 n 2 


+ 


C 2 ?li 

16n 2 


4c 2 nf 
96 n 2 


\c\ \c\ri 2 c 

8 niri 2 96 n\ 256nf' 


c 2 n,2 


C 2 7^2 


Xu = 7 r n 2 


1 + T-j-(n - 3n 2 ) + -2 (4 r 


—24nri2 + SOrj.^)] , 


Xi,2 = 8 tt 2 ?ti 


4 4 

1 + -j —-An — 2m) + — (2?r 2 


+ 10 ?i 2 — 12 nni)] . 


(D3) 


384n 2 


(D4) 


(D5) 


(D 6 ) 


In the above calculation for compressibility, the condition 
H = 2(m — p 2 ) + c 2 /4 was used, i.e. dH = 0 gives (up 
to the O (^ 2 ) order): 


valid in ID interacting systems because there does not 
exist a well-defined quasiparticle. Here we demonstrate 
that exact solution of the TBA equations does show a 
novel existence of Fermi liquid like signature in ID in¬ 
teracting systems. In fact, the additivity rules which we 
found in previous sections reflect a ID Fermi liquid like 
nature, also see a discussion in Oil- 

Feedback interaction. The excitations near Fermi 
points in ID many-body systems can form a collective 
motion independent of microscopic details of systems, 
namely, there exists a certain dispersion relation between 
energy and momentum. We observe that integrable sys¬ 
tems provide a deep understanding of the intrinsic con¬ 
nection between the TLL and the Fermi liquid. In fact, 
the TBA equations of ID systems give the exact dressed 
energies and determine the dispersion relations of each 
branch. At low temperatures, only the behavior of the 
dressed energies near the zero point, or saying, the fci? 
and vf at the ID Fermi points, determines the first or¬ 
der and second orders of thermodynamic quantities. In 
general, for the attractive SU{N ) Fermi gases, the spin 
fluctuations are suppressed. Therefore the TBA equa¬ 
tions can be rewritten as 


dni _ U 2 _ f _ 8 n 2 6 ni n\ _ 16 n\ \ 

dn 2 4ni \ |c| |c| 2 |c|ni |c|n 2 y ’ 

Moreover, the interaction effect enters into the collective 
velocities v\, V 2 of the excess fermions and bound pairs. 
For strong attraction, they are given by m 

Vi S3 -^27rm (l + 8 n 2 /|c| + 48n\/c 2 ) , 

2 m 

V 2 ~ -^- 7 m 2 (l + 2A/|c| + 3A 2 /c 2 ) , (D 8 ) 

2 m 

with A = 2iii + n 2 . 


Appendix E: Connection between the Luttinger and 
Fermi liquids 


Wilson Ratios. Now we further build up a connec¬ 
tion between the TLL and the Fermi liquid. By definition 
([I]) and ([2]), the two type of Wilson ratios of interacting 
Fermi liquid in 3D are given by 


px _ _ p« _ _ 

"w- l+Fir Hw ~ 1 




(El) 


which depends on Landau parameters Fq’ s charactering 
the interaction ; 3|. This is very similar to our finding 
for ID systems, see Eqs. (12) and (13). However, it’s 
extremely hard to calculate the Landau parameters F, q’“ 
in Fermi liquid theory due to the reason that N* (0) can 
not be obtained explicitly. Fermi liquid theory elegantly 
maps an interacting system into a free fermion system 
where the interaction is encoded into the density of state 
and effective mass. But the Fermi liquid theory is not 


N 

e {r) (k) =e°(k) - ^ A rs * (r = 1,2,- • • ,w), 

S= 1 

min(m,n ) 

A-mn = ^ ^ Q"m+n—q (E2) 

q =1 


where 


^0 _ h- 


2 m r{k 2 - &.) = 2rm 


— P n 




is the first 


order coefficient describing excitation energy of a single 
r-complex, here /i r the effective chemical potential. Here 
Po = rhk. The function a m (x) is defined in Appendix 
A. These equations show a similar form of ’feedback in¬ 
teraction’ equation in the Landau Fermi liquid theory 
rm rn. This encourage us to further find a perturba¬ 
tion idea into our calculations. What below gives a Fermi 
liquid like description for the low energies of the ID at¬ 
tractive Fermi gases, i.e. mapping an interacting system 
with multi-subsystems to a multicomponent free system. 

Phenomenological description. The symmetry 
group of the interaction of quasi particles in the ID 
“Fermi liquid” is reduced from the SO( 3) group to the 
cyclic group C 2 . Therefore, interaction parameter be¬ 
tween the ’quasi’ momentum p and p' could be written 
as = /g’ a + fi’ a sign(p-p/). Then following the con¬ 
ventional Fermi liquid theory, we still can have Landau 
parameters Fq which are interacting parameters [16J 
and describe the main properties of our system. Conven¬ 
tional sound velocity is associated with oscillations in the 
density of a fluid, and hydrodynamics gives that 


P 


K 

mn 


(E3) 


Where p = mn is the density of the fluid and k = — 
is the bulk modulus. In FL theory n = (1 + Fq), 
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where the IV* (0) denotes the density of state in momen¬ 
tum space near Fermi surface, therefore the velocity of 
first sound is given by 


v 


2 


n 

mN* (0) 


(l + *o) 


(E4) 


By the ID analog of Fermi liquid theory, we have the 
relations [H] n = IV* (0) = and m = m*( 1 + 
F-f). Then we obtain: 

v 2 3 4 5 6 7 8 9 10 =v 2 F (l + F$)(l + Ff). (E5) 


ample, the velocities of trions, pairs and single fermionic 
atoms, are given by 


3 1 (3) 2 1 (2) 1 (1) 

- — 'll' '• 'll - - 'll' '* 'll - 'll' ' 


v = 3 « f 


W = 2 Vp 


V = Vr 


(E7) 


respectively. This builds up an intrinsic connection be¬ 
tween the Luttinger liquid and Fermi liquid. 

Furthermore, from the free fermion nature of the Fermi 
liquid, we can express the specific heat in terms of these 
sound velocities 


Here, due to the collective motion, the backflow does not 
exist in ID, therefore Ff « 0. This can be seen from the 
dressed energy equation (E2), where the effective mass 
rm is almost unchanged up to the order of 0 (c~ 2 ) for a 
strong coupling regime. 

Consistency. Using the exact solutions, we can cal¬ 
culate the Wilson ratio R ^ for different phases of TLLs. 
For a single state of a r-complex, the Wilson ratio is 
known from Eq. (El I. Then the parameter Fq x could be 
determined by a comparison with the Bethe ansatz result 
in strong coupling limit (i.e. |c| —> oo) via 


R 


c _ 

w — 


1 

¥+Ff) 


{ 1 ; for free fermions 

4; for pairs . (E6) 

9; for trions 



7T / 1 1 1 

3 \td v 2 + v 3 


,(E8) 


where we have dropped the unit , here the interac¬ 

tion effect is encoded in the velocities v 1,2,3 . This equa¬ 
tion could be obtained by consiering the leading term 
of the TBA equations at low temperature region. This 
shows a consistency of our derivation by the ID Fermi 
liquid like description. The derivation above could be ex¬ 
tended to the SU(w) attractive Fermi gases in a straight¬ 
forward way. All the derivation above could be directly 
extended to the SU(w) attractive Fermi gases: 


Cy 7T / 1 1 1 

— = -- 1 - 1 - 1 - 

T 3 yu 1 v 2 v u 


(E9) 


From the equation (E51 and (E6) we could calculate the 
relation between Fermi speed and sound velocity. For ex- 


We can also prove this additivity rule using the TBA 
equations, see Appendix B. 
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